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ABSTRACT 

Cold fronts in cluster cool cores should be erased on short timescales by thermal conduction, unless 
protected by magnetic fields that are "draped" parallel to the front surfaces, suppressing conduction 
perpendicular to the sloshing fronts. We present a series of MHD simulations of cold front formation 
in the core of a galaxy cluster with anisotropic thermal conduction, exploring a parameter space of 
conduction strengths parallel and perpendicular to the field lines. Including conduction has a strong 
effect on the temperature distribution of the core and the appearance of the cold fronts. Though 
magnetic field lines are draping parallel to the front surfaces, preventing conduction directly across 
them, the temperature jumps across the fronts are nevertheless reduced. The geometry of the field is 
such that the cold gas below the front surfaces can be connected to hotter regions outside via field 
lines along directions perpendicular to the plane of the sloshing motions and along sections of the front 
which are not perfectly draped. This results in the heating of this gas below the front on a timescale 
of a Gyr, but the sharpness of the density and temperature jumps may be nevertheless preserved. By 
modifying the gas density distribution below the front, conduction may indirectly aid in suppressing 
Kelvin-Helmholtz instabilities. If conduction along the field lines is unsuppressed, we find that the 
characteristic sharp jumps seen in Chandra observations of cold front clusters do not form; therefore, 
the presence of of cold fronts in hot clusters is in contradiction with our simulations with full Spitzer 
conduction. This suggests that the presence of cold fronts in hot clusters could be used to place upper 
limits on conduction in the hulk of the ICM. Finally, the combination of sloshing and anisotropic 
thermal conduction can result in a larger flux of heat to the core than either process in isolation. 
While still not sufficient to prevent a cooling catastrophe in the very central (r ~ 5 kpc) regions of 
the cool core (where something else is required, such as AGN feedback), it reduces significantly the 
mass of gas that experiences a cooling catastrophe outside those small radii. 

Subject headings: conduction — cooling flows — galaxies: clusters: general — X-rays: galaxies: 
clusters — methods: hydrodynamic simulations 



1. INTRODUCTION 

X-ray observations of galaxy clusters often show a 
bright, dense, centrally peaked core comprised of the in- 
tracluster medium (ICM) much cooler than the cluster 
average. The X-ray radiative cooling time, inversely pro- 
portional to the square of the gas density, is very short in 
these dense cores, on the order of several hundred Myr in 
the very central regions (jSarazinlll986l^ . This gave rise to 
an early "cooling flow" scenario, in which the gas cools 
catastrophically in the cluster center and is continuously 
replaced by the inflow of the surrounding cooling gas. 
However, Chandra and XMM- Newton have failed to de- 
tect the implied large quantities of very cool gas in the 
cores; in particular, the characteristic emission lines such 
as Fe XVI I emitted by the plasma w ith T < 1 keV are not 
observed (|Peterson fc Fabianll2QQ6l ). Since X-ray cooling 
is directly observed, a mechanism is required to heat the 
ICM in the core to compensate for the cooling. The two 
most studied heating mechanisms are feedback from ac- 
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tive galactic nuclei (AGN) and thermal conduction. In 
this paper, we will be discussing the latter mechanism in 
connection with the phenomenon of gas sloshing. 

The role of thermal conduction in the context of heat- 
ing the cool cores of clusters of galaxies has been stud- 
ied by many authors ( see, e .g., 'Binnev fc Cowid *19811;^ 
Naravan fc MedvedevI I2QQ1 : [Ruszkowski fc Begelman 
2QQ2I : IZakamska fc NaravanI [2003. : .Guo et al.. ^ ,2008. ). 



Since in the ICM the Larmor radius of the electrons spi- 
raling along the B ~ 0.1 — 1 fiG magnetic field lines 
is many orders of magnitude smaller than the electron 
mean free path, the conduction is strongly anisotropic, 
occurring essentially exclusively along the field lines. 
Therefore, most of these works assumed an isotropically 
tangled magnetic field geometry and paramaterized the 
heat flux as an effective isotropic thermal condu ctivity 
given as a fraction, /sp, of the ideal (field- free) iSpitzerl 
(|1962r ) heat fiux. 

However, theoretical and numerical work over the last 
decade has revealed that the stability and thermal prop- 
erties of the ICM are possibly more complicated than 
previously thought if anisotropic thermal conduction is 
m odelled explicitly in the MHD regime. Early analysis 
bv lBalbusI (| 2 0001 . 1 2 00 If ) showed that when magnetized at- 
mospheres with anisotropic thermal conduction are con- 
sidered, their stability depends on the temperature gradi- 
ent rather than the entropy gradient. He demonstrated 
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that in regions where the temperature increases in the 
direction of gravity, the plasma is susceptible to the mag- 
netothermal instability (hereafte r MTI) . 
Further analytical work by iQuataertI (|2QQ8[) and 

numerical wo rk by iParrish fc Q uataert (2008) and 

[Parrish etUH (120081, 120091) has shown that the cluster 
atmosphere is unstable even when the temperature de- 
creases in the direction of gravity, as in the cluster cool 
cores (heat-flux buoyancy instability; HBI). The effect 
of the HBI is to realign the magnetic fields azimuthally 
(6 • f 0), and thereby significantly reduce the effective 
radial heat flux Qr- This renders thermal conduction 
almost useless against preventing a "cooling catastro- 
phe." Subsequent work has shown that relatively low 
levels of turbulence can rearrange the field lines in a 
more isotropic c onfiguration, permitting heat conduc- 
tion t o the core (jParrish et alJ l2Q10il : iRuszkowski fc Otil 



I2OIOI . I2OII). More recent l y, and of direct relevance to 
this work, iLecoanet et al.l (|2012f ) simulated the dynam- 



ics of Ray leigh- Taylor (RT) stable and unstable contact 
discontinuities with anisotropic thermal conduction. Of 
particular relevance is their finding that the tempera- 
ture gradients of RT-stable contact discontinuities (such 
as sloshing cold fronts) can be smoothed by anisotropic 
thermal conduction if the field lines are at least partially 
misaligned with the surface of the discontinuity. The 
smoothing of the temperature gradient of the RT-stable 
contact discontinuities in these simulations was limited 
by the HBI, which reoriented magnetic field lines ini- 
tially perpendicular to the front surface in a more par- 
allel direction, suppressing conduction across the surface 
at later times. 

So far, the numerical simulations of the evolution of 
the cluster atmospheres with anisotropic thermal con- 
duction have mostly been limited to isolated, spherically 
symmetric clusters initially in hydrostatic equilibrium. 
These works have been essential in order to determine 
the character of the instabilities that are present when 
such physical effects are taken into account. A recent ex- 
ception to this is the paper by Ruszkowski et al. (2011), 
which simulated the formation of a galaxy cluster from 
cosmo logical initial conditions. 

The previous works assumed that the heat conduction 
along the field lines is at the Spitzer value. However, the 
effective conductivity along the field lines may be reduced 
due to strong curvature of the field line s and stochastic 
magnetic mirrors on small scales (e.g ., iChandran et al. l 
0,999; Malvshkin fc Kulsrud 2001; Nar avan fc M edvedev 
[2001) . Markevitc h et al.i (^2003b) attempted to constrain 
heat conduction in the ICM from the observed tem- 
perature differences in the merging, non-cool-core clus- 
ter Abell 754. They determined that the effective con- 
duction in the direction of the observed temperature 
gradients is an order of magnitude less than Spitzer, 
which, under the assumption of an isotropically tan- 
gled field, imples suppres s ion a long the field lines. Fi- 
nally, Schekochihi n et aTl (j2005l ) suggested that plasma 
microscale instabilities may be very important in the 
high-/3 plasma of the ICM, and consequently its thermal 
conductivity may be m any orders of magnitude smaller 
than the Spitzer value (Schek ochihin et al]l2008f ). 

In this work we consider a scenario where the ini- 
tial cluster core is relaxed but undergoes a period of 



gas "sloshing." Many cool-core systems exhibit edges 
in X-ray surface brightness approximately concentric 
with respect to brightness peak of t he cluster (e.g., 
iMazzotta et"alll200ll : iMarkevitch efaIll200lL [20031 ). X- 
ray spectra of these regions have revealed that in most 
cases the brighter (and therefore denser) side of the edge 
is the colder side, and hence these jumps in gas den- 
sity hav e been dubbed "cold fronts" ( for a detailed re- 
view see IMarkevitch fc VikhlininI 120071 ) . These features 
have been observed in the majority of cool-core systems 
(Ghizzardi et al. 2010). The presence of one or more cold 
fronts is seen as an indication of the motion of gas in a 
cluster. For relaxed systems, the primarily spiral-shaped 
cold fronts are believed to arise from gas sloshing in the 
deep dark matter-dominated potential well. These mo- 
tions are initiated when the low-entropy, cool gas of the 
core is displaced from the bottom of the dark matter po- 
tential well, either b y gravitational per turbations from 
infalling subclusters (jAscasibar fc Mark evitch 2006) or 
by a n interaction with a shock front (jChurazov et al.l 

In an earlier paper, ZuHone et ajj (|2010l . hereafter 
ZMJIO), we have investigated the possible flow of heat to 
the cluster cool core as a result of the sloshing motions 
bringing hot gas from larger radii into contact with the 
cool gas of the core. We determined that if the gas viscos- 
ity is small, significant mixing between the hot and cold 
gas will result, raising the entropy of the core, though 
this heat flow was found to be insufficient to completely 
offset the effects of radiative cooling. In that work, we 
did not include thermal conduction, but only mixing of 
gases with different temperatures. 

ZMJIO modeled the ICM as a (possibly viscous) un- 
magnetized fluid. However, the ICM is known to be mag- 
netized, and this fact has important considerations for 
heat transport. For example, it has been suggested that 
the narrow width of cold fronts in galaxy clusters can 
be used to constrain the effectiveness of heat conduc- 
tion i n the ICM (e.g..'Ettori fc Fabian '20001 'Xiang et al.l 
I2OO7I ). The sharp temperature gradient of the cold front 
should be smoothed out on short timescales if conduction 
is not suppressed across the front surface. Such suppres- 
sion could be provided by magnetic field lines oriented 
parallel to the cold front surfaces, a situation believed 
to naturally arise due to the effects of m agnetic draping 
and shear amphfication ('Asai et al."2004. 120051 : ILvutikovl 
^06; Asai et al. 2007; Keshet et al. 20101). In a previous 
paper (?, hereafter ZMLll), we performed simulations 
of gas sloshing with magnetic fields, showing that such 
shear amplification and draping does occur around cold 
fronts in cluster cores. These simulations also showed 
that magnetic fields have similar effects as viscosity on 
the ICM, inhibiting the mixing of the hot and cold flows 
and reducing the increase of core entropy due to sloshing. 

In this work, we build upon our earlier simulations 
(ZMJIO, ZMLll) by including the effects of anistropic 
thermal conduction in a simulation of gas sloshing in a 
cool-core galaxy cluster. We will show that the magnetic 
flelds draped across the cold fronts indeed inhibit the 
smoothing out of the temperature gradient, but only to 
a degree. The crucial points to note are a) due to fluid 
instabilities the magnetic fleld lines may not be oriented 
perfectly azimuthal with the front and b) cold fronts are 
not completely enclosed surfaces-there are still fleld lines 
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that connect the cold side of the cold fronts to regions of 
higher temperature in the cluster atmosphere over short 
distances. Therefore, if thermal conduction is important, 
it may provide a more efficient way than gas mixing to 
transfer heat between the hot and cold flows character- 
istic of sloshing motions. 

This paper is organized as follows: in Section 2 we de- 
scribe the simulations and the code. In Section 3 we de- 
scribe the effects of anisotropic conduction and gas slosh- 
ing on the thermal state of the cluster core for different 
configurations of the magnetic field and various prescrip- 
tions for thermal conduction. Finally, in Section 4 we 
summarize our results and discuss future developments 
of this work. We assume a flat ACDM cosmology with 
h 0.7, 0.3, and 0.02/^-^. 

2. METHOD 

In our simulations, we solve the following set of MHD 
equations (in Gaussian units): 



dp 

m 



V • (pv)=0 



dt 
dE ^ 



pvv - — ) + Vp = pg 



(1) 
(2) 
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where pth is the gas pressure, e is the gas internal energy 
per unit volume, and B is the magnetic field strength. 
For all our simulations, we assume an ideal equation 
of state with 7 = 5/3 and primordial abundances with 
molecular weight /i = 0.6. 

We performed our simulations using FLASH 3.3, a 
parallel hydrodynamics/ A/'-body astrophysical simula- 
tion code developed at the FLASH Center for Computa- 
tiona l Science at the Uni versity of Chicago (iFrvxell et al.l 
120001 : iDubev et al"]l2009[ ). FLASH 3.3 solves the equa- 
tions of magnet ohydro dynamics using a directionally 
unsplit staggered mesh algorithm (USM; iLee fc Deand 
i009). The USM algorithm used in FLASH 3.3 is based 
on a finite- volume, high-order Godunov scheme com- 
bined with a constrained transport method (CT), which 
guarantees that the evolve d magnetic f ield satisfies the 
divergence- free condition (jEvans fc Hawlev ^988). In 
our simulations, the order of the USM algorithm cor- 
responds to the P iecewise-Parabolic Method (PPM) of 
IColella fc Woodw ard (1984), which is ideally suited for 
capturing shocks and contact discontinuties (such as the 
cold fronts that appear in our simulations). 

The heat flux due to thermal conduction is given as: 



Q 



/^(l-66)+/||66 



■ VT 



(7) 



where we chose to represent the conductivity coefficient 
as a product of the Spitzer conductivity n<^^ and the frac- 
tions /|| and of Spitzer conductivity parallel and per- 
pendicular to the field line, respectively, hh and 1— hh 
are projection operators that single out the components 
of the temperature gradient VT parallel and perpendic- 
ular to the magnetic field line, with h the unit vector 
in the direction of the local magnetic field. Most works 
assume f\\ = 1 and f±=Om the conditions of the ICM 
(e.g., heat conducted strictly along the field lines without 
suppression), due to the large mean free path for colli- 
sions compared to the Larmor radius of the electron. We 
assume this prescription as the default, but also inves- 
tigate alternative scenarios, represented by simulations 
with different values of /y and f±. For example, if the 
magnetic field is strongly curved on small scales, the ef- 
fective conductivity along the field lines may be smaller 
than the Spitzer value, resulting in a smaller /m. 

Alternatively, if a fraction of the magnetic field lines 
draped along the cold front surfaces reconnect with lines 
below these surfaces, a small amount of heat conduction 
perpendicular to the cold front surfaces may result. Our 
ideal MHD simulations do not model reconnection ex- 
plicitly, as non-ideal terms including resistivity are not 
included (there is a small numerical resistivity associated 
with the finite resolution of the simulation, but this effec- 
tive resistivity is negligible, see Section ??). However, we 
may account for the conduction that may occur across 
cold front surfaces in an approximate way by modeling it 
as a non-zero f± , with the same temperature dependence 
as the parallel conduction. Since we are not modeling the 
reconnecting field lines and the conduction along them 
explicitly, our simple model is intended to provide only 
an order of magnitude estimate of this possible effect. 

It is important to distinguish our "perpendicular" con- 
duction across the cold front surfaces from conduction 
perpendicular to the field lines resulting from collisions 
between electrons. Because the Larmor radius is many 
orders of magnitude smaller than the electron mean free 
path, this perpendicular conduction is effectively zero 
for the conditions in galaxy clusters. Additionally, the 
temperature dependence of this conduction would be 
(X T~^l^ instead of the ( x T^^'^ dependen ce for conduc- 
tion along the field lines (|Braginskiilll965f ). In contrast, 
since our "perpendicular" conduction is intended to ac- 
count for conduction that occurs parallel to reconnecting 
field lines, it has the same temperature dependence as the 
parallel con duction. 

Following ICowie fc McKed (|1977[ ). we included the ef- 
fect of conduction saturation whenever the characteris- 
tic lengthscale associated with the temperature gradi- 
ent is smaller than the mean free path, though in the 
bulk of the ICM this effect is not significant. We imple- 
mented anisotropic c onduction following the approach of 
iSharma fc Hammettl (j2007). More specifically, we ap- 
plied a monotonized central (MC) limiter to the conduc- 
tive fiuxes. This method ensures that anisotropic con- 
duction does not lead to negative temperatures in the 
presence of steep temperature gradients. A verification 
test of our implementation of the a nisotropic conduc- 
tion m odule in FLASH 3.3 was given in lRuszkowski etld] 
(|20TT1 ). 

We want to examine the effects of anisotropic ther- 
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mal conduction in a sloshing galaxy cluster core in iso- 
lation (in particular, the effect on the cold fronts), and 
so for most of our simulations, we do not include the 
effects of radiative cooling. We then include radia- 
tive cooling to determine the effect of sloshing with 
anisotropic conduction on the thermal state of the cluster 
core. The X-ray cooling function is determined using the 

ewe. Kaastra. fc Liedahll ()1995f ) model with a uniform 
metallicity of Z = O.SZq. 

The gravitational potential on the grid is set up as the 
sum of two "rigid bodies" corresponding to the contribu- 
tions to the potential from both clusters. This approach 
to modeling the potential is used for simplicity and speed 
over solving the Poisson equation for the matter distribu- 
tion, and is an adequate approximation for our purposes. 
It is the same approach that we used in ZMLll, and is 
justified in Roediger & ZuHon^ (j2Qlir ). 

We refer the reader to ZMLll for the full details of 
our physical setup, but will provide the following short 
description. Our simulations consist of a massive (M ^ 
1.5 X lO^^A^o, Tx ^ S keV), cool-core cluster, initially 
in hydrostatic equilibrium, merging with a small (mass 
ratio R = 5) gasless subcluster, which sets off the sloshing 
of the cool core. We begin the simulation at a point in 
time when the cluster centers have a mutual separation 
of d = 3 Mpc and an impact parameter b = 500 kpc. The 
initial velocities of the subclusters are set up assuming 
that the total kinetic energy of the system is set to half of 
its total potential energy. For all of the simulations, we 
set up the main cluster within a cubical computational 
domain of width L = 2.4 Mpc on a side, with a finest 
cell size on our AMR grid of Ax = 2.34 kpc. 

The magnetic field of the cluster is set up in an 
identical way to t he sim ulation s of ZMLll, fo llow- 
ing the approach of iRuszko wski et al.l ()2007l l2008f ) and 
iRuszkowski fc OE (|2010l ). A random magnetic field B(k) 
is set up in k-space on a uniform grid using indepen- 
dent normal random deviates for the real and imaginary 
components of the field. Thus, the components of the 
complex magnetic field in k-space are set up such that 

B,{k) = Bi[N{ui)^iN{u2)] (8) 

By{k) = B2[N{us)^iN{u^)] (9) 

B,{k) = Bs[N{u,)^iN{ue)] (10) 

where N{u) is a function of the uniformly distributed 
random variable u that returns Gaussian-distributed ran- 
dom values, and the values Bi are field amplitudes. We 
adopt a dependence of the magnetic field amplitude B{k) 
on the wavenumb e r |k| simil ar to (but not the same a s) 
lR.uszkowski et al.l (j2007[ ) and lR.uszkowski fc OhI (j2010f ): 

B{k) (X k-^^^^exp[-{k/kof]exp[-ki/k] (11) 

where ko and ki control the exponential cutoff terms 
in the magnetic energy spectra. The cutoff at high 
wavenumber ko roughly corresponds to the coherence 
length of the magnetic field /cq = 27r/Ao = 27r/(43 kpc), 
observed in some cluster cores. The cutoff at low 
wavenumber ki = 27r/Ai roughly corresponds to Ai ~ 
^50o/2 ~ 500 kpc, which is held fixed for all of our simula- 
tions. This field spectrum corresponds to a Kolmogorov 
shape for the energy spectrum (Pe(/c) ex k~^^^) with 
cutoffs at large and small linear scales. This field is then 



Fourier transformed to yield B(x), which is rescaled to 
have an average value of y^87rp/(3 to yield a field with a 
pressure that is proportional to the gas pressure, i.e., to 
have a spatially uniform P for the initial field. We set 
the initial P = 400. This procedure ensures that there 
will be significant power at all length scales between Ai 
and Ao, which will result in a field that is significantly 
tangled down to the scale of Aq. 

As mentioned above, previous simulation works have 
predicted that due to the action of the HBI, an azimuthal 
magnetic field should arise in the cluster core, in the ab- 
sence of turbulence driven by galaxy motions or AGN. 
For this purpose, we also include a simulation that has 
a purely azimut hal initial field similar to that used in 
Bogda novic et al.l (j2009i) . described in spherical coordi- 
nates by 

Br = (12) 
5^ = 2^0 sin l9 cos 20 (13) 
B^ = -Bo sin 2(/) sin 20 (14) 

Its magnitude was then scaled by the gas pressure with 
P = 400. 

Finally, for both configurations, a field consistent with 
V • B = is generated by performing the Fourier trans- 
form and then by "cleaning" the field of divergence terms 
in Fourier space via the the projection 

B'(k) = (I-kk)B(k) (15) 

where k is the unit vector in k-space, and I is the identity 
operator. Note that this operation does not change the 
power spectrum of the magnetic field fluctuations. Upon 
transformation back to real space, the magnetic field sat- 
isfies V • B = 0. This field is then interpolated onto our 
AMR grid in such a way that the condition V • B = is 
maintained. 

Because our sloshing simulations begin with the two 
clusters at a distance, it takes a period of time (~2 Gyr) 
for the influence of the subcluster to perturb the cool core 
of the main cluster and initiate sloshing. If we had be- 
gan the simulation with anisotropic thermal conduction 
"switched on", it would alter the thermal state of the 
core and the HBI and MTI would alter the cluster mag- 
netic field spatial configuration long before sloshing had 
begun, since the growth times of these instabilities and 
the conduction timescale for the core are on the order 
of ~ a few 10^ yr. Worse yet, if we began the simu- 
lation with radiative cooling switched on, the cool core 
would cool rapidly (also on a timescale of ~ a few 10^ yr) 
and we would have a cooling catastrophe before sloshing 
started. To avoid these effects, we switch on conduction 
(and cooling in the simulations where it is included) at 
a time t ^ 2 Gyr after the beginning of the simulation, 
shortly after the subcluster has made its pericentric pas- 
sage of the cluster core and right at the beginning of the 
sloshing period. This enables us to examine the effects of 
conduction and sloshing together on a typical cool-core 
cluster state. 

An unfortunate side-effect of this evolutionary period 
before the core passage is that the magnetic field config- 
uration will slowly evolve away from its initial configura- 
tion. Hence, the initially isotropically tangled field con- 
figuration with (1 6 • f|) = 0.5 evolves to a slightly more 
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TABLE 1 
Simulation Parameter Space 



Simulation 


Relaxed / Sloshing 


Magnetic Field 


f.. 

h 


J± 


X-ray Cooling? 


SI 


Sloshing 


Tangled 


0.0 


0.0 


NO 




Sloshing 


Azimuthal 


U.U 


U.U 


IN 


SCI 


Sloshing 


Tangled 


1.0 


0.0 


NO 


SC2 


Sloshing 


Azimuthal 


1.0 


0.0 


NO 


SC3 


Sloshing 


Tangled 


0.1 


0.0 


NO 


SC4 


Sloshing 


Tangled 


0.1 


0.01 


NO 


RX 


Relaxed 


Tangled 


0.0 


0.0 


YES 


RCX 


Relaxed 


Tangled 


1.0 


0.0 


YES 


sx 


Sloshing 


Tangled 


0.0 


0.0 


YES 


SCXl 


Sloshing 


Tangled 


1.0 


0.0 


YES 


SCX2 


Sloshing 


Azimuthal 


1.0 


0.0 


YES 


SXNoB 


Sloshing 


No Fields 


0.0 


0.0 


YES 



t = 2.25 Gyr SI 


S2 

1 


SCI 


~ SC2- 

"J 


SC3 


SC4 



9.0 



7.5 



> 

c r\ ^ 
6.0 ^ 



4.5 



'3.0 



Fig. 1. — Slices through the center of the domain in the x — y plane of the temperature for the sloshing simulations with and without 
conduction for the epoch t = 2.25 Gyr. Each panel is 500 kpc on a side. The colorscale is temperature in keV. Major tick marks indicate 
100 kpc distances. 



radial field configuration with {\b • f|) ^ 0.56, and the 
initially azimuthal field configuration with {\b • f|) = 
also evolves to a slightly more radial field configuration 
with {\b • f|) ~ 0.35. This is not a major concern, as the 
former is not too far off from an isotropically tangled field 
configuration and the latter is not too far from the field 
configuratio n of the saturation state of the HBI as seen i n 
simulations ( Parrish et al.|[2QQ9l : lBogdanovic et al.ll2009l ). 
In real clusters, there is always turbulence and small 
disturbances from galaxi es and DM clumps that keep 
the field tangled (Parrish et al.|[2QTQl : fRuszkowski & Oh 
l20iai2QllD. 

We will refer to the simulations using the following 
codes. The letter "R" will indicate a cluster that is 
evolved in isolation, whereas "S" will indicate a cluster 
evolved with a gravitational interaction causing slosh- 



ing. "C" will indicate anisotropic thermal conduction is 
turned on, and "X" will indicate X-ray radiative cooling 
is turned on. Simulations that differ in the prescription 
for conduction and the initial field line configuration will 
be delimited by numbers. Table [1] lists all simulation 
runs with their particular parameters. 

3. RESULTS 

3.1. Brief Outline of the Sloshing Process 

We have described the initiation and the evolution 
of sloshing cold fronts in previous works (e.g., AM06, 
ZMJIO, ZMLll), the latter paper including the effects 
of magnetic fields. We refer the reader to these papers 
for the full description, but we will outline the process in 
brief here. 

As the subcluster approaches the main cluster's core 
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Fig. 2. — Same as Figure [T] but corresponding to the epoch t = 2.75 Gyr. 



and makes its passage, the gas and DM peaks of the main 
cluster feel the same gravity force toward the subcluster 
and move together towards it. However, the gas feels 
the effect of the ram pressure of the ambient medium. 
This fact becomes significant as the gas core is held back 
from the core of the dark matter by this pressure. This 
results in a separation of the gas and DM peaks of the 
cluster. As this ram pressure declines with the decreas- 
ing gravitational attraction of the outgoing subcluster, 
the gas that was held back by the ram pressure falls into 
the DM potential minimum and overshoots it. In addi- 
tion to the gravitational disturbance, the wake trailing 
the subcluster transfers some of the angular momentum 
from the subcluster to the core gas and also helps to 
push the core gas out of the DM potential well. As the 
cool gas from the core climbs out of the potential min- 
imum, it expands adiabatically. However, the densest, 
lowest-entropy gas quickly begins to sink back towards 
the potential minimum against the ram pressure from 
the surrounding ICM. Once again, as the cool gas falls 
into the potential well it overshoots it, and the process 
repeats itself on a smaller linear scale. Each time, a con- 
tact discontinuity ( "cold front" ) is produced. Due to the 
angular momentum transferred from the subcluster by 
the wake, these fronts have a spiral-shaped structure. As 
we noted in ZMLll, this general outline is not affected 
significantly by the strength of the magnetic field or its 
spatial configuration, so the features seen in density and 
temperature maps are essentially the same between our 
two reference (e.g., without conduction) simulations, SI 
and S2 (the SI simulation corresponds to the BetaJ^OO 
simulation in ZMLll), and the runs without magnetic 
fields in ZMJIO. When speaking of the orientation of 
the temperature gradient across the cold fronts, we shall 
refer to the hotter, less dense gas as the side "above" 
the front and the colder, denser gas "below" the front. 



consistent with the direction of the gravitational acceler- 
ation. 

3.2. The Effect of Conduction on Cold Fronts 

Figures [T] through [H show shces through the center of 
the cluster of the gas temperature for several different 
epochs, for our simulations with and without anisotropic 
conduction (but without radiative cooling). We will 
first compare runs with no conduction to runs with full 
Spitzer conduction along the field lines (/y = l,/± = 0). 
Our runs SI (top left panels) and S2 (top center panels) 
correspond to sloshing simulations with initially tangled 
fields and initially azimuthal fields, respectively. Like- 
wise, our runs SCI (top right panels) and SC2 (bot- 
tom left panels) correspond to sloshing simulations with 
anisotropic heat conduction, with initially tangled fields 
and initially azimuthal fields. 

Including conduction has a dramatic effect on the tem- 
perature structure of the cold fronts. The adiabatic simu- 
lations without conduction {SI and S2) develop a number 
of well-defined cold fronts over time, which maintain tem- 
perature differences of ~ 1.5-2 across the fronts through- 
out their evolution. This is not the case for the simula- 
tions with anisotropic unsuppressed Spitzer conduction 
(simulations SCI and SC2). Over time, the cold gas in- 
side the fronts is heated up, so that the temperature con- 
trasts across the fronts become smaller, roughly ~ 1.3- 1.5. 
The changing thermal state of the gas also changes the 
shape of the fronts in some places. In particular, due to 
the fact that the inner 50 kpc becomes nearly isothermal, 
cold fronts do not occur on the smallest scales, in con- 
fiict with observations of cool-core clusters, which typ- 
ically show a number of cold fronts at both small and 
intermediate radii within the core. The temperature of 
the core in simulation SC2 increases more slowly than in 
simulation SCl^ because the more azimuthally oriented 
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field lines in the former suppress conduction from larger 
radii more efficiently than the more tangled fields of sim- 
ulation SCI do. 

Figures [T] through |4] also show results for the simu- 
lations with other values for the conduction coefficients 
and /||. The SC3 simulation has the same field line 
configuration as the SCI simulation, but the conduction 
along the lines is suppressed by a factor of 10. If field 
lines are strongly curved on very small linear scales, the 
effective conduction coefficient parallel to the field lines 
may be smaller than the Spitzer value by an order of mag- 
nit udedueto^t he presence of stochastic magnetic mirrors 
fe.g.. iNaravan Medvedevll2QQlf ). We choose /y = 0.1 
to determine the effect of such a reduction. Alterna- 
tively, if plasma microscale instabilities play a significant 
role, the thermal conductivity of the ICM may be man y 
orders of magnitude smaller ( S cheko chihin et a l? ' 2 8 V 
For all practical purposes, this case is equivalent to the 
runs without conduction. 

On the other hand, if (for example) ~1% of the mag- 
netic field lines parallel to the fronts reconnect, the ef- 
fective fraction of Spitzer conduction perpendicular to 
the large-scale direction of the field could be f±_ ^ 0.01. 
We assume this value of f± for simulation 5(7^, keep- 
ing the value of /y = 0.1 for comparison with simulation 
SC3. As expected, for the SC3 simulation, far less heat 
is conducted to the cluster core, and the temperature 
contrast of the cold fronts with the surrounding medium 
is much closer to the case without any conduction. In 
the SC4 simulation, the small perpendicular conductiv- 
ity, combined with a steep temperature jump across the 
cold front, has succeeded over time to smooth out the 
temperature gradient across the cold fronts. 

Figure [5] compares temperature gradients across cold 
fronts between the SI and the SCI simulations at a late 
epoch, t = 3.25 Gyr. Four particular locations across the 
fronts are picked out for inspection (see Figure [3] for the 
profile locations, marked by the corresponding letters). 
The profiles show that the temperature jumps in the 
conducting simulation have all been reduced. For some 
of the profiles (e.g., (a) and possibly (d)), the width of 
the temperature jump in the simulation with conduction 
SCI is not greater than the jump in the non-conducting 
simulation SI (a few simulation cells, ~ 5 — 10 kpc- 
essentially unresolved by the simulations), whereas for 
other profiles in simulation SCI, the temperature jumps 
have been smoothed out over several zone sizes (we will 
explore these issues further in Section 13. 3p . Particularly 
apparent in profile (d) is that conduction has modified 
the temperature distribution below the front, while pre- 
serving the sharp jump across the front. 

Figure [5] also highlights example temperature jumps 
in the SC3 and SC4 simulations (along the same pro- 
file locations, see Figure [3]). The temperature difference 
across the fronts are similar between these two simula- 
tions and the SI simulation, due to the high suppression 
of conduction, but the smoothing out of the gradient 
due to the conduction perpendicular to the field lines in 
the SC4 simulation is noticeable, with the temperature 
change occurring over roughly twice the number of cells 
(~10 kpc) compared to the case without perpendicular 
conduction due to field line reconnection. 

The jumps in density across the cold fronts are simi- 
larly affected by conduction. As heat flows from the hot 



gas above the cold fronts to the cool gas below the front, 
the former contracts and the latter expands in the same 
proportion as the change in temperature. Figure [6] shows 
the density jumps across the cold fronts at the same po- 
sitions in Figure [H In simulations with conduction, the 
density jumps are reduced and the widths of the density 
jumps are increased to varying degrees, in the same pro- 
portion as the temperature jumps except that the sign 
of the gradient is reversed. In the most extreme cases, 
the density jumps have been essentially erased, either in 
the strong parallel conduction case SCI or the case SC4 
with a small component of conduction perpendicular to 
the magnetic field lines. The smoothed profiles in these 
simulations are in disagreement with observations of cold 
fronts in clusters. 

The differences between the simulations SI, SC3, and 
SC4 indicate that numerical reconnection across the cold 
front surfaces and any associated heat flux due to this 
reconnection is negligible. The cold fronts in the simula- 
tions SI and SC3 have essentially the same temperature 
and density jumps, whereas adding a small perpendicular 
conduction smooths out the jumps over twice the num- 
ber of cells. The fact that we do not see this degree of 
smoothing in the SC3 simulation indicates that numer- 
ical reconnection is too small to provide any significant 
heat flux across the fronts. 

The smaller density jumps and the modification of the 
density profiles below the fronts that reduce the density 
contrast have a significant effect on the appearance of 
cold fronts as seen in X-ray images, due to the strong 
dependence of the X-ray emission on the gas density. 
Figure [3 shows synthetic X-ray images for the SI, SCI, 
and 5(7^ simulations, at the epoch t = 2.75 Gyr, created 
by computing the bolometric MeKaL emissivity in each 
FLASH cell and projecting it along the z-axis (in the 
same manner as ZMLll). Additionally, a Chandra X- 
ray image of Abell 2319 is also presented for comparison 
(this cluster is slightly hotter than our model cluster, 
with average T ~ 9 — 10 keV). This cluster exhibits a cold 
front with a very sharp jump in X-ray emi ssion, with a 
width roughly the size of the Chandra PSF (jO'Hara et al.l 
2004: Govoni et al 2004). 

The characteristic spiral shape of the cold fronts is 
apparent in each simulated image. However, the sharp 
edges in emission across the front surfaces are essentially 
absent from the SCI case, whereas in the SC3 and SI 
simulations, the sharp jumps are clearly seen. There- 
fore, for our simulated hot cluster, unsuppressed Spitzer 
conduction appears to produce cold fronts that are far 
less prominent than those in observed hot clusters, and 
it is questionable whether or not they may be even called 
"cold fronts" at all. We reserve a more quantitative com- 
parison with observed clusters of different temperatures 
and plasma /3 values and the possible constraints on con- 
ductivity for a future paper. 

It also is apparent from these figures that conduction 
may have an effect on the smoothness of the front sur- 
faces. For example, the cold fronts in simulations SCI 
and SC2 exhibit fewer small-scale features along the 
front surfaces than their counterparts without conduc- 
tion, SI and S2. The widening of the cold front inter- 
faces (seen clearly in Figure [6]) will cause them to be less 
susceptible to the Kelvin-Helmholtz instability (see, e.g., 
.Churazov fc Inogamov.2004) . 
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Fig. 5. — Temperature profiles of cold fronts in simulations without conduction and with varying prescriptions for conduction, along the 
profiles marked with the corresponding letters in Figure [3] Conduction reduces the magnitude and increases the width of the jumps to 
varying degrees. 
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3.3. The Geometry of the Magnetic Field Lines 

The simulations from ZMLll showed magnetic field 
is amplified and the lines oriented parallel to the cold 
front surfaces, which would prevent conduction across 
the cold front surfaces from the hotter gas just above 
the fronts. The fact that there are examples of discern- 
able temperature jumps over the same spatial range (a 
few simulation cells) even in the SCI run indicates that 
conduction is still strongly suppressed directly across the 
front surfaces. What accounts for the significant changes 
in density and temperature across the fronts that still 
occur? 

The cool gas below the cold fronts is only completely 
immune to the effects of heat conduction if there are no 
field lines leading from hotter regions of the cluster at- 
mosphere into the cold regions of the fronts from other 
directions. The fronts seen in X-ray observations do not 
completely enclose the cool core. In principle, there could 
be magnetic field lines that reach the regions below the 
fronts from hotter regions on the other side of the cluster. 



More importantly, simulations of gas sloshing show that 
the resulting spiral structure has an alternating pattern 
of hot and cold fiows, with hot flows de veloping below 
cold fronts as they expand to higher radii. iKeshetl (|2Qllf ) 
showed that these latter flows are of hot gas being drawn 
into the core region form higher radii. These regions can 
be magnetically connected to the cool gas immediately 
below the fronts. Additionally, though the cold fronts 
that form in simulations are extended in directions per- 
pendicular to the sloshing plane, the characteristic size of 
the fronts in this direction (along the line of sight where 
the observer sees the spiral structure) is smaller than in 
the plane. This exposes the cold regions of the fronts to 
magnetic fleld lines extending into hotter regions along 
directions perpendicular to the sloshing plane. 

The preceding arguments lead us to expect that some 
fleld lines that pass through cold gas below the front sur- 
faces may be connected to hotter regions over relatively 
short distances. Figures [9] and [10] show magnetic fleld 
lines beginning at several points in the cold region below 
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Fig. 7. — Projected X-ray emission along the z-axis of the simulation domain for the SI (no conduction), SCI (Spitzer conduction), and 
SC3 (0.1 Spitzer conduction) simulations at the epoch t = 2.75 Gyr, with a Chandra X-ray image of A2319 included for comparison. White 
bars indicate 100 kpc distances. This indicates that conduction smooths out the density jumps, making them barely discernable as seen in 
X-ray images in comparison to the sharp jumps in emission seen in A2319. 



the front surface and followed for a distance of 100 kpc. 
A large number of field lines are tangentially oriented to 
the front surfaces, as can be easily seen in both projec- 
tions, shielding this cold gas from the hotter gas directly 
above the front. However, the figures show that there 
are still a number of field lines that extend into hotter 
regions of the cluster core along other directions. At 
the epochs of t = 2.5 Gyr and t = 3.0 Gyr, the coldest 
temperatures below the front surfaces are approximately 
T ^4.5-5 keV, and these regions are connected to regions 
with temperatures of T ~6-7.5 keV. In the right panel 
of Figure [9l field lines stretched along the y and z direc- 
tions connect cold gas underneath the front surfaces to 
hotter gas. In Figure [TOl lines connecting cool gas with 
hot gas that was brought in by sloshing are evident near 
the center of the left panel of the figure. 



Finally, the simulations show that there is not always 
a perfect alignment of strong magnetic field layers with 
cold fronts. Figure [8] shows slices through the center 
of the simulation domain in the x — y plane (coincident 
with the plane of the sloshing motions) of the gas tem- 
perature and the magnetic field energy for the epochs t = 
2.75 Gyr and t = 3.25 Gyr with the magnetic field vector 
components in the x and ^/-directions overlaid. Extended 
regions of strong magnetic field indicate places where the 
magnetic field strength has been increased by shear am- 
plification. The locations of the profiles from Figure [3] 
are also marked. While some locations along the cold 
fronts exhibit a strong field layer parallel to the front, 
this is not strictly true everywhere along the fronts. In 
particular, the t = 2.75 Gyr panel shows that the trail- 
ing edge of the cold spiral in the southwestern direction 
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Fig. 8. — Slices through the z = coordinate plane of the temperature and the magnetic field energy in the vicinity of four representative 
profiles with the 2D vector field b = (bx, by) overlaid for the simulation SCI (Spitzer conduction along field lines) at the epochs t = 2.75 Gyr 
and 3.25 Gyr. Each panel is 400 kpc on a side. Tick marks indicate 100 kpc distances. The field vectors are located at the positions of 
their tails. Top panels: Temperature in keV. Bottom panels: The magnetic field energy B'^ in /xG^. Red, elongated regions in the bottom 
panels indicate magnetic draping layers. The white rectangles indicate the regions without strong draping layers described in the text. The 
white lines indicate the profile locations from Figure [3] 



(indicated roughly by the white rectangles in Figure [8j) 
is not draped by a high-5 layer. It can be seen that 
the magnetic field vectors along this portion of the spiral 
front are not strictly aligned parallel to the front surface. 
Regions of ordered, strong field are interspersed with re- 
gions of weak, tangled field as shown by the vectors at 
these locations. At a somewhat later time, t = 3.25 Gyr, 
this location still lacks a strongly magnetized parallel 
layer. This also happens to be the location of profile 
(d). That particular temperature jump was shown to 
be smoothed out significantly in Figure [5l Though at 
the epoch t = 3.25 Gyr there does appear to be a small, 
parallel, highly magnetized layer at that particular front 



location, this layer was not present previously and the 
temperature jump likely has been smoothed out before 
the formation of this layer. 

Cold fronts that are not draped by strong magnetic 
fields will be more susceptible to the K-H instability, 
which will distort the front surface and lead to gas mix- 
ing between hot and cold phases. This is because the 
orientation of the draping fields tends to coincide with 
the direction of the shearing motions. Perturbations of 
the field lines perpendicular to the front surfaces at these 
locatio ns will lead to sm oothing out of the fronts as out- 
lined in lLecoanet et al.l ((2012), though the extent of this 
smoothing will be limited by the HBI, which will re- 
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Fig. 9. — Magnetic field lines in the SCI (Spitzer conduction) simulation beginning at points inside the cold fronts with T < 4.5 keV and 
continuing up to a length of 100 kpc at the epoch t = 2.5 Gyr. Left panel: View along the 2;-axis. Right panel: View along the x-axis. 
Color along the lines indicates temperature in keV. The box size is 150 kpc, and the major tick marks indicate 50 kpc distances. This 
indicates that field lines connect regions of cold gas under the cold front surfaces to hot gas outside them. 
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Fig. 10. — Magnetic field lines in the SCI simulation beginning at points inside the cold fronts with T < 5.5 keV and continuing up to a 
length of 100 kpc at the epoch t = 3.0 Gyr. Left panel: View along the z-axis. Right panel: View along the x-axis. Color along the lines 
indicates temperature in keV. The box size is 250 kpc, and the major tick marks indicate 50 kpc distances. This indicates that field lines 
connect regions of cold gas under the cold front surfaces to hot gas outside them. 



orient magnetic field lines parallel to the front surface. 
iLecoanet et al.l (|2012) estimated that the scale height of 
the temperature gradient for such a front would smooth 
to a height of Ht ^ y^^^ /q^^^^ where x is the Spitzer 
thermal diffusivity and g is the gravitational accelera- 
tion. Taking example values from our simulation for a 
cold front at r ^120 kpc, we find x 5 x 10^^ cm^ s~^ 
and ^ ~ 3 X 10~^ cm s~^. This yields an approximate 
scale height of the temperature gradient across the front 
of Ht ~ 30 kpc, which is similar to the widths of some of 
the temperature gradients across the fronts in Figure [5l 
In simulations of thermal instabilities in magnetized 
galaxy clusters, a helpful quantity to examine has been 



the average angle between the magnetic field and radial 
unit vectors, given by 



9b) =cos-i(|6-f|) 



(16) 



Isotropically oriented tangled field lines will have an av- 
erage (1 6 • f|) = 0.5 and hence the corresponding angle 
is {9b) ^ 60°. The HBI drives field lines in a more az- 
imuthal direction {Ob ~ 70 — 80°) in the presence of a 
positive radial temperature gradient (as in a cool core). 
Sloshing motions in general preserve the overall positive 
temperature gradient of the cool core, implying that the 
HBI should still be active. However, they will also drag 
the field lines around into a more azimuthally oriented 
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Fig. 11. — The evolution of the volume-averaged (within a sphere 
of radius R = 100 kpc) angle between the magnetic field unit vec- 
tor and the radial unit vector {Ob) = cos~-'^{|6 • r\) over time for 
simulations with and without conduction. 



configuration as tiiey are draped around tiie cold fronts. 

Figure [11] shows the evolution of the volume- averaged 
angle of the magnetic field with respect to the radial di- 
rection within a radius of 100 kpc over time for the Sl^ 
S2^ SCl^ and 5(7^ simulations. The simulations with ini- 
tially tangled and azimuthal fields begin with {Ob) ^ 56° 
and {Ob) ^ 70°, respectively. We find that largely the 
evolution of this "magnetic angle" is driven by the slosh- 
ing motions themselves, since the evolution of 6>b is essen- 
tially identical between the simulations with and with- 
out conduction that have the same initial magnetic field 
setup. An exception to this general rule are the simula- 
tions with initially azimuthal magnetic field lines, S2 and 
SC2. When conduction is included, the final Ob indicates 
a more azimuthal direction than in the case where con- 
duction is absent. This is possibly due to the fact that 
in this case the core is more magnetically isolated due 
to the azimuthally oriented field lines, and the positive 
temperature gradient is better preserved in the core than 
in the tangled field line case. In such a scenario, the HBI 
would drive the field lines in a more azimuthal direction 
once the conduction is turned on. 

3.4. The Evolution of the Entropy of the Cluster Core 

Given that our simulations show that the magnetic 
field geometry of the sloshing cool core still permits heat 
conduction to the cluster core, it is important to deter- 
mine what effect this has on the thermal state of the 
core. The evolution of average entropy within the inner 
50 kpc of the cluster center over time is shown in Figure 
[T2I In the simulations without conduction {SI and S2)^ 
the average entropy of the cluster core increases due to 
sloshing bringing hot gas from the cluster outskirts with 
the cool gas of the core into contact, which then mix to 
a certain degree. However, the magnetic field suppresses 
such mixing (ZMLll), and the increase in entropy levels 
off roughly 0.5 Gyr after the simulation begins. ZMLll 
showed that the difference in initial field geometry has 
little effect on the entropy increase in this case. 
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Fig. 12. — Evolution of the average entropy within 50 kpc of the 
cluster center for simulations with and without anisotropic thermal 
conduction. 

In contrast, in the simulations with anisotropic Spitzer 
conduction {SCI and SC2)^ heat conducts along the 
field lines to the cluster center, and the entropy increase 
of the core continues unabated for the duration of the 
simulation. As expected, the effect is stronger in the 
case with initially random field lines than in the case 
with initially azimuthal field lines. The average spe- 
cific entropy increase in these two simulations is roughly 
~ 20 — 30 keV cm^ higher than in the adiabatic 
simulations. 

A similar situation exists for the simulations with 
weaker conduction coefficients. The simulation with con- 
ductivity 1/10 of the Spitzer value strictly along the field 
lines {SC3) has a correspondingly smaller increase in en- 
tropy over the SCI simulation, with ~ 5 keV cm^, 
close to the entropy increase due to sloshing alone, but 
the entropy is still increasing at the end of the simulation. 
The inclusion of even the small, 1/100 Spitzer conduc- 
tion perpendicular to the field lines {SC4) results in a 
significantly higher entropy increase over simulation SI 
than the SC3 case, with AS' ~ 10 keV cm^-this is be- 
cause in the sloshing core, temperature gradients across 
the magnetic field lines are much higher than along the 
lines. The result of all of the simulations with some form 
of anisotropic conduction is that heat can be conducted 
efficiently to the core, even in the presence of azimuthal 
cold fronts. 

3.5. Simulations with Conduction and Radiative 
Cooling 

We now include radiative cooling in our simulations. 
The runs with radiative cooling SCXl and SCX2 corre- 
spond to the non-cooling runs SCI and SC2. The cooling 
is "turned on" at the same moment as the anisotropic 
conduction, at the epoch t = 2 Gyr. As control cases 
to isolate the effects of sloshing and conduction, we also 
ran the following simulations. RX is an initially relaxed 
cluster left in isolation without heat conduction. Due 
to the absence of any heating mechanism, it experiences 
a cooling catastrophe in the core within a few hundred 
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Fig. 13. — Slices through the center of the domain in the x — y plane of the temperature for the sloshing simulations with cooling for the 
epoch t = 2.25, 2.75, 3.25, and 3.75 Gyr. Each panel is 500 kpc on a side. The colorscale is temperature in keV. Major tick marks indicate 
100 kpc distances. 
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Myr. RCX is a similar run with an isolated cluster, ex- 
cept it also includes Spitzer conduction along the field 
lines. Thermal conduction partially offsets the cooling of 
the core in this run, but the HBI quickly reorients the 
magnetic field lines in the core in the azimuthal direc- 
tion and a cooling catastrophe still occurs. The last run, 
SX, is identical to simulation SI (without conduction and 
with sloshing), but with radiative cooling included. 

In ZMJIO, we showed that radiative cooling may be 
slowed by sloshing itself, due to the mixing of hot gas 
from the outskirts with the cold gas for the core. How- 
ever, the inclusion of even weak magnetic fields reduces 
this mixing and inhibits the consequent heating of the 
core, depending on the strength of the magnetic field 
(ZMLll). From our results in ZMJIO, we expect that 
regardless of the strength of the magnetic field, the heat 
delivered to the core from sloshing alone will not be suf- 
ficient to prevent a cooling catastrophe at the innermost 
regions of the cool core. The left panels of Figure [T3l 
shows sHces of temperature through the center of the 
cluster for selected epochs for the SX simulation. In this 
run, the core temperature drops below 1 keV within a ra- 
dius of r ~ 20 kpc. Additionally, the core density peaks 
at high values of ng ~ 1 — 5 cm~^, in confiict with obser- 
vations of real clusters. Therefore, the goal of the runs 
with conduction and cooling is to determine if sloshing 
with conduction together is able to prevent such a cooling 
catastrophe. 

The center and right panels of Figure [13] shows slices 
of temperature through the center of the cluster for se- 
lected epochs for the sloshing simulations including cool- 
ing and conduction. Anisotropic conduction, even in the 
case where it is most efficient, is unable on its own to 
suppress the runaway cooling in the very center of the 
cluster, with central temperatures falling to the temper- 
ature ffoor of our simulation (T ~ 0.01 keV), far below 
the temperatures observed in real clusters. However, the 
temperature is significantly raised outside the very cen- 
tral region of the cool core, at radii r ^ 10 kpc. 

The two leftmost panels of Figure [TH show the evolu- 
tion of the gas mass with short cooling times within the 
inner 50 kpc of the cluster center over time. We were 
only able to run the RX simulation out to an epoch of 
t ^ 3.6 Gyr, as the extremely high densities, low temper- 
atures, and magnetic pressures in this simulation made 
the simulation too unstable to continue past this point. 
Finally, we have also included a cooling run SXNoB^ 
identical to simulation SX, but without magnetic fields, 
to compare the rate of increase of cool gas when the 
mixing of hot and cool gases is not suppressed by such 
fields. This the same simulation as in ZMJIO, but the 
setup (resolution and model for the gravitational poten- 
tial) is identical to our present work, to enable a direct 
comparison. 

The top panels of Figure [TH show the evolution of the 
gas mass with short cooling times within the inner 50 kpc 
of the cluster center over time. We were only able to 
run the RX simulation out to an epoch oi t ^ 3.6 Gyr, 
as the extremely high densities, low temperatures, and 
magnetic pressures in this simulation made the simula- 
tion too unstable to continue past this point. Finally, 
we have also included a cooling run SXNoB, identical to 
simulation SX, but without magnetic fields, to compare 
the rate of increase of cool gas when the mixing of hot 



and cool gases is not suppressed by such fields. This the 
same simulation as in ZMJIO, but the setup (resolution 
and model for the gravitational potential) is identical to 
our present work, to enable a direct comparison. 

Over the course of ~2 Gyr, all of the simulations have 
a significant increase of the mass of cool gas in their 
cores. However, both sloshing and conduction reduce 
the buildup of cool gas. The top-left panel of Figure [T4l 
shows the buildup of gas mass with tcooi < 0.1 Gyr. At 
the onset of cooling, there is no gas in the core with a 
cooling time this short. Within ~0.25 Gyr, in the re- 
laxed cluster simulations {RX and RCX)^ a gas mass of 
~ 10^ — 10^^ Mq with such cooling time has built up 
in the core. This is somewhat slowed by conduction in 
the RCX simulation, and by t ~ 4 Gyr the mass of cool 
gas is lower than that of the RX simulation by a fac- 
tor of a few. In the simulations with sloshing [SX and 
SCXl)^ the initial buildup of cool gas does not begin un- 
til ~0.5 Gyr after the onset of cooling. However, in the 
SX simulation, the mass of cool gas eventually catches 
up with that of the RX simulation. Out of all of these 
simulations, the run with both sloshing and conduction 
(SCXl) is most effective at suppressing the cooling of 
the gas for the long term, with nearly an order of mag- 
nitude less mass of cool gas in the core than in the RX 
simulation. 

In the case of simulation SXNoB, no gas below tcooi = 
0.1 Gyr appears until ~1 Gyr after the onset of cooling. 
Gas mixing, which is much more efficient in the absence 
of a magnetic field (or perhaps in the presence of field 
variations on microscopic scales that would strongly re- 
duce field tension), provides a source of heat to the core 
gas and prevents the catastrophic cooling of gas until 
this later stage. However, the mass of cool gas increases 
quickly after this point and is eventually comparable with 
that of simulation RCX. This indicates that in the ab- 
sence of magnetic field tension, strong gas mixing due 
to sloshing is a comparable source of heat to anisotropic 
conduction. 

The top-right panel of Figure [14] shows the buildup of 
gas mass with tcooi < 1-0 Gyr. In this case, all sim- 
ulations already start with the same mass of gas with 
cooling timescales this short. However, the general trend 
for each of these simulations is the same as for the case 
for the gas mass with tcooi < 0.1 Gyr. The runs with 
conduction (or gas mixing) and sloshing exhibit a slower 
buildup of cool gas. 

The bottom panels of Figure [14] show slices through the 
z = plane of the temperature at the epoch t = 3.5 Gyr, 
with contours marking the surfaces with specific cooling 
times of 0.1, 0.5, and 1.0 Gyr, for the simulations RX 
and SCXl. It is clear that the effect of sloshing and con- 
duction together is to reduce the volume of the region 
with very low cooling times and temperatures. With- 
out conduction and sloshing, a r ~ 10 kpc region with 
^cooi < 0.1 Gyr develops, whereas with both of these ef- 
fects included, this region is restricted to r ~ 3 kpc. 

Overall, this result indicates that sloshing, even with 
anisotropic thermal conduction, is not enough to com- 
pletely stabilize a cool core against a cooling catastrophe, 
and that for the central regions of clusters, AGN feedback 
or other mechanisms are necessary to regulate the tem- 
perature of the core. However, sloshing and conduction 
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Fig. 14. — Top panels: Mass of gas with a cooling time below a given value vs. simulation time for several of the cooling simulations. 
Left panel: Mass of gas with tcool < 0.1 Gyr. Right panel: Mass of gas with tcool < 1-0 Gyr. Bottom panels: Slices through the z = 
plane of the temperature in keV for the RX and SCXl simulations at the epoch t = 3.5 Gyr. White contours mark the surfaces with tcool 
= 1.0, 0.5, and 0.1 Gyr. Each panel is 75 kpc on a side. Tick marks indicate 10 kpc distances. 



are able to suppress or significantly delay catastrophic 
cooling in the cluster cooling cores (r 50 — 100 kpc) 
outside the very central regions of radii r ^ 3 kpc. 

4. DISCUSSION AND CONCLUSIONS 

Gas sloshing, a prevalent phenomenon in relaxed 
galaxy clusters, is evidenced by the nearly ubiquitous 
presence of spiral-shaped cold fronts in cluster cool cores. 
The temperature jumps across these fronts imply that 
either thermal conduction is intrinsically weak in galaxy 
clusters or the magnetic field is oriented parallel to the 
front surfaces, restricting conduction across the fronts. 
Most previous works have assumed the latter, and sim- 
ulations have confirmed that such magnetic layers form 
around cold fronts in clusters. 

Sloshing may also have an important role in determin- 
ing the thermal evolution of the cluster core. It brings 
the cold gas of the core into contact with hot gas, possi- 
bly facilitating a transfer of heat between these phases, 
either via mixing of the hot and cold gas or by heat con- 
duction. This raises the possibility that sloshing may be 
partially responsible for preventing a "cooling catastro- 
phe" in the cores of galaxy clusters. However, the ICM 



is magnetized, which, as we showed in ZMLll, partially 
suppresses mixing. 

To determine whether anisotropic heat conduction in 
the presence of magnetic fields has any effect on sloshing 
cold fronts and the thermal state of the cool cores, we 
have performed MHD simulations with various configu- 
rations of the magnetic field and the conduction coeffi- 
cients. We determine that, in spite of the fact that the 
cold fronts indeed are mostly draped by a strong par- 
allel magnetic field, they are not completely thermally 
isolated from hotter gas. The shear amplification and 
stretching of magnetic field lines occurs only in a thin 
layer at some of the cold front surfaces, which only pre- 
vents heat conduction to the cold gas below the front 
from the hot gas directly above the front. Heat can still 
fiow around the cold front, along field lines connected 
to the cold gas from other regions, either from the hot 
gas at higher radii perpendicular to the sloshing plane 
or from the hot fiows that are brought inside the cool 
core by sloshing. Thus, while the sharpness of the tem- 
perature gradient at the cold front surfaces may be pre- 
served at some locations along the fronts, the cool gas 
below the fronts can still be heated from other direc- 
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tions. Also, at some points along the cold front surfaces, 
the sharpness of the temperature gradient is not pre- 
served because a stabilizing parallel magnetic field layer 
fails to form. Conduction also leads to the reduction of 
the density jumps and the widening of the cold front in- 
terfaces, which can aid in suppressing Kelvin-Helmholtz 
instabilities. Importantly, for the hot simulated clus- 
ter presented in this work, we find that if conduction 
is unsuppressed along the field lines, the density jumps 
can be smoothed to such a degree that the cold fronts 
are no longer apparent in synthetic X-ray observations, 
in profound disagreement with observed clusters. Due 
to the strong temperature dependence of Spitzer con- 
duction, we cannot use our model cluster, which is hot 
(Tx ~ 8 keV), to draw conclusions about the smearing 
out of cold fronts due to conduction in cooler clusters. In 
a future paper, we will attempt more quantitative com- 
parisons with observed cold fronts in other clusters in 
order to more generally constrain conduction. 

We find that a conduction coefficient parallel to the 
field lines an order of magnitude smaller than the Spitzer 
value has a nearly insignificant effect on the temperature 
of the core and the density and temperature jumps at 
the cold front surfaces, producing cold fronts that are in 
agreement with observations. On the other hand, adding 
even a small perpendicular conduction (0.01 of Spitzer), 
that may, for example, mimic the effect of field line recon- 
nection, has an effect on the core specific entropy similar 
to the case with unsuppressed Spitzer conduction along 
the field lines (because the temperature gradients across 
the field lines are higher), and slightly smoothes out the 
temperature gradients of the fronts over the course of a 
few Gyr. 

In the absence of radiative cooling, gas mixing and the 
conduction of heat between the hot and cold phases in 
the cluster core raises the temperature and the entropy 
of the core, and if left unchecked, would eliminate the 
cool core entirely, leaving behind a high-entropy isother- 
mal core. When radiative cooling is taken into account, 
we find that, even for full Spitzer conduction along the 
field lines, a cooling catastrophe still occurs in the clus- 
ter center, though the effects of conduction and sloshing 
together can significantly suppress the buildup of cool 
gas in the core and restrict the cooling catastrophe to a 
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small (r ^ 5 kpc) volume in the very center. The fact that 
conduction is unable to suppress cooling completely in 
the cluster center is somewhat inevitable, as for isobaric 
perturbations the cooling rate increases with lower tem- 
perature and the heating rate from conduction decreases 
(jKunz et al. 2011). AGN feedback or other mechanisms 
are still necessary to prevent gas from cooling in these 
central regions. However, our results show that slosh- 
ing and conduction together can produce a significant 
amount of heat to offset some cooling in the cluster core. 

We note that we have only considered one merger sce- 
nario and one initial P for the magnetic field. Though 
different choices in these parameter spaces are not ex- 
pected to change our main conclusions, from our previous 
works we may make reasonable guesses about the effect 
of changing our initial conditions. Certain merger setups 
result in more mixing of hot and cold gases than others 
(ZMJIO), whereas a lower initial P results in stronger 
magnetic field layers at the cold front surfaces and hence 
more suppression of K-H instabilities (ZMLll). Most im- 
portantly for our considerations, the strong temperature 
dependence of Spitzer conduction makes it particularly 
strong for our T ~ 8 keV simulated cluster. In order to 
use observations of cold fronts in real clusters to put con- 
straints on suppression of conduction, a wider parameter 
space of initial cluster temperatures and plasma p values 
would need to be explored. Such a study, including a 
quantitative comparison with a sample of observed cold 
fronts, will be presented in a future work. 



JAZ thanks Ian Parrish and Mikhail Medvedev for use- 
ful discussions and advice. Calculations were performed 
using the computational resources of the National In- 
stitute for Computational Sciences at the University of 
Tennessee and the Advanced Supercomputing Division 
at NASA/ Ames Research Center. Analysis of the simu- 
lation data was carried out using the AM R analysis and 
visualization toolset yt (|Turk et al.ll2011l ). which is avail- 
able for download at http://yt-project.org. JAZ is 
supported by the NASA Postdoctoral Program. MR ac- 
knowledges NSF grant 1008454. The software used in 
this work was in part developed by the DOE NNSA-ASC 
OASCR Flash Center at the University of Chicago. 



Cowie, L. L., & McKee, C. F. 1977, ApJ, 211, 135 
Dubey, A., Antypas, K., Ganapathy, M. K., Reid, L. B., Riley, 
K. M., Sheeler, D., Siegel, A., Weide, K. Extensible component 
based architecture for FLASH, a massively parallel, 
multiphysics simulation code. Parallel Computing 35 (10-11), 
512-522. 

Ettori, S., & Fabian, A. C. 2000, MNRAS, 317, L57 
Evans, C. R., & Hawley, J. F. 1988, ApJ, 332, 659 
Fryxell, B., et al. 2000, ApJS, 131, 273 

Ghizzardi, S., Rossetti, M., & Molendi, S. 2010, A&A, 516, A32 
Govoni, F., Markevitch, M., Vikhhnin, A., et al. 2004, ApJ, 605, 
695 

Guo, F., Oh, S. p., & Ruszkowski, M. 2008, ApJ, 688, 859 
Keshet, U., Markevitch, M., Birnboim, Y., & Loeb, A. 2010, ApJ, 

719, L74 
Keshet, U. 2011, arXiv: 1111.2337 

Kunz, M. W., Schekochihin, A. A., Cowley, S. C, Binney, J. J., & 

Sanders, J. S. 2011, MNRAS, 410, 2446 
Lecoanet, D., Parrish, I. J., & Quataert, E. 2012. larXiv:1202. 21171 
Lee, D., &; Deane, A. E. 2009, Journal of Computational Physics, 

228, 952 



19 



Lyutikov, M. 2006, MNRAS, 373, 73 
Malyshkin, L., &; Kulsrud, R. 2001, ApJ, 549, 402 
Markevitch, M., Vikhlinin, A., & Mazzotta, P. 2001, ApJ, 562, 
L153 

Markevitch, M., Vikhlinin, A., & Forman, W. R. 2003, 

Astronomical Society of the Pacific Conference Series, 301, 37 

Markevitch, M., Mazzotta, P., Vikhlinin, A., et al. 2003, ApJ, 
586, L19 

Markevitch, M., & Vikhlinin, A. 2007, Phys. Rep., 443, 1 
Mazzotta, P., Markevitch, M., Vikhlinin, A., Forman, W. R., 

David, L. P., & VanSpeybroeck, L. 2001, ApJ, 555, 205 
Mewe, R., Kaastra, J. S., & Liedahl, D. A. 1995, Legacy, 6, 16 
Narayan, R., & Medvedev, M. V. 2001, ApJ, 562, L129 
O'Hara, T. B., Mohr, J. J., & Guerrero, M. A. 2004, ApJ, 604, 

604 

Parrish, I. J., & Quataert, E. 2008, ApJ, 677, L9 
Parrish, I. J., Stone, J. M., & Lemaster, N. 2008, ApJ, 688, 905 
Parrish, I. J., Quataert, E., &; Sharma, P. 2009, ApJ, 703, 96 
Parrish, I. J., Quataert, E., &; Sharma, P. 2010, ApJ, 712, L194 
Parrish, I. J., McCourt, M., Quataert, E., & Sharma, P. 2012, 
MNRAS, 2543 

Peterson, J. R., & Fabian, A. C. 2006, Phys. Rep., 427, 1 
Quataert, E. 2008, ApJ, 673, 758 
Roediger, E., &; Zuhone, J. A. 2011, MNRAS, 1793 
Ruszkowski, M., &; Begelman, M. C. 2002, ApJ, 581, 223 
Ruszkowski, M., Enfilin, T. A., Briiggen, M., Heinz, S., & 
Pfrommer, C. 2007, MNRAS, 378, 662 



Ruszkowski, M., EnBhn, T. A., Briiggen, M., Begelman, M. C, & 

Churazov, E. 2008, MNRAS, 383, 1359 
Ruszkowski, M., & Oh, S. P. 2010, ApJ, 713, 1332 
Ruszkowski, M., & Oh, S. P. 2011, MNRAS, 414, 1493 
Ruszkowski, M., Lee, D., Briiggen, M., Parrish, I., & Oh, S. P. 

2011, ApJ, 740, 81 
Sarazin, C. L. 1986, Reviews of Modern Physics, 58, 1 
Schekochihin, A. A., Cowley, S. C, Kulsrud, R. M., Hammett, 

G. W., & Sharma, P. 2005, ApJ, 629, 139 
Schekochihin, A. A., Cowley, S. C, Kulsrud, R. M., Rosin, M. S., 

& Heinemann, T. 2008, Physical Review Letters, 100, 081301 
Sharma, P., & Hammett, G. W. 2007, Journal of Computational 

Physics, 227, 123 
Spitzer, L. 1962, Physics of Fully Ionized Gases, New York: 

Interscience (2nd edition), 1962 
Turk, M. J., Smith, B. D., Oishi, J. S., Skory, S., Skillman, S. W., 

Abel, T., &; Norman, M. L. 2011, ApJS, 192, 9 
Vikhhnin, A., Markevitch, M., & Murray, S. S. 2001, ApJ, 549, 

L47 

Vikhlinin, A. A., & Markevitch, M. L. 2002, Astronomy Letters, 
28, 495 

Xiang, F., Churazov, E., Dolag, K., Springel, V., &; Vikhlinin, A. 

2007, MNRAS, 379, 1325 
Zakamska, N. L., & Narayan, R. 2003, ApJ, 582, 162 
ZuHone, J. A., Markevitch, M., & Johnson, R. E. 2010, ApJ, 717, 

908 (ZMJIO) 

ZuHone, J. A., Markevitch, M., & Lee, D. 2011, ApJ, 743, 16 



